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Data from solar neutrino and KamLAND experiments have led to a discovery of nonzero neutrino 
masses. Here we investigate what these data can tell us about neutrino interactions with matter, 
including the poorly constrained flavor-changing i/e — Vt interactions. We give examples of the 
interaction parameters that are excluded by the solar/KamLAND data and are beyond the reach of 
other experiments. We also demonstrate that flavor-changing interactions, at the allowed level, may 
profoundly modify the conversion probability for neutrinos of energy < 6 MeV and the values of the 
mass parameter inferred from the data. The implications for future experiments are discussed. 



I. INTRODUCTION 

For several decades, experiments have been trying to 
test the Standard Model (SM) paradigm that neutrinos 
are massless and interact only via the W and Z gauge 
boson exchange. In recent years, a breakthrough has 
emerged: data from the solar, atmospheric, and reac- 
tor neutrino experiments have indicated that neutrinos 
do have masses, and hence the SM is incomplete. It is 
the right time to ask whether the SM predictions for the 
neutrino-matter interactions can be similarly tested. 

The aim of this Letter is to investigate what can be 
learned about neutrino-matter interactions from present 
and future solar and KamLAND neutrino data. We an- 
swer two questions: (i) Can the solar and KamLAND 
experiments constrain parts of the parameter space that 
are presently inaccessible by non-oscillation experiments? 
(ii) Can the uncertainty in our present knowledge of 
neutrino-matter interactions affect the determination of 
the oscillation parameters? As we show, the answer 
to both questions is affirmative. We give explicit ex- 
amples of parameters that are disfavored by solar and 
KamLAND data and that are beyond the reach of non- 
oscillation experiments. We also demonstrate that non- 
standard interactions (NSI) , at an allowed level, can qual- 
itatively modify the fit to the data and change the values 
of inferred mass parameters. This scenario leads to non- 
trivial predictions for future experiments. A full presen- 
tation of the numerical constraints we obtain is beyond 
the scope of this Letter and will be given elsewhere 0. 

We study the effects of NSI on top of nonzero neu- 
trino masses and mixing, since the latter are required by 
the KamLAND data. This scenario has been previously 
studied "5, '3, 1^ with an emphasis on the limit of small 
flavor-changing interactions. We extend the formalism 
developed in these papers to make it applicable to our 
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problem. 



II. NSI AND SOLAR NEUTRINOS: THE 
PHYSICS 

Low-energy neutrino interactions can be described by 
four-fermion interaction vertices, L 3 ^ipipip. The 
vertices affecting neutrino evolution in matter are those 
containing two neutrino lines (2i/). In the SM, these 
vertices receive contributions from neutral current (NC) 
processes and, if the initial state contains a charged lep- 
ton, also charged current (CC) processes. The NC pro- 
cesses are predicted to be flavor-preserving and universal. 
Possible non-standard (both flavor-preserving and flavor- 
changing) contributions to the 2^ vertices can, most gen- 
erally, be parameterized as 

L^S/ ^ -2V2Gj.(P„7pi^/3)(e^^'^/L7''/L+e^r/fl7''//?) 
+ h.c. (1) 

Here e^^'^ {^{//' ) denotes the strength of the NSI between 
the neutrinos v of flavors a and /3 and the left-handed 
(right-handed) components of the fermions / and f^Gp 
is the Fermi constant. 

Bounds on the epsilons come from accelerator-based 
experiments, such as NuTeV |6| and CHARM 0, and 
experiments involving charged leptons. In the later case, 
we do not include bounds obtained by the SU{2) sym- 
metry, since strictly speaking these can be avoided if, 
for example, the corresponding operators contain Higgs 
doublets 1^. Both types of experiments are quite effec- 
tive at constraining the vertices involving the muon neu- 
trino, giving eef,< 10"^, e^^< 10"^ - 10"^ At the same 
time, bounds on See, Cer, and err are rather loose, e.q., 
|e™^| < 3, -0.4 < e™^ < 0.7, |e™| < 0.5, |e^^| < 0.50. 

NSI can modify both the neutrino propagation (oscil- 
lation) [1^] and neutrino detection processes. The prop- 
agation effects of NSI are, first of all, only sensitive to 

e{/p when f = f ^ (henceforth, e^Jf = e^^), and, sec- 
ond, only to the vector component of that interaction. 
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= ^ap + ^Q^- The matter piece of the osciUation 
Hamihonian can be written (up to an irrelevant constant) 
as: 
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where rig is the number density of electrons in the 
medium. The epsilons here are the sum of the contri- 
butions from electrons (e*^), up quarks (e"), and down 
quarks (e'') in matter: = Z]/=n.d,e ^a/3"'//"e- Hence, 
unlike in the standard case [tap = 0), the NSI matter ef- 
fects depend on the chemical composition of the medium. 

The CC detection reactions at SNO, KamLAND and 
the radiochemical experiments, just like the production 
reactions in the Sun, are unchanged by Eq. JQ). On the 
other hand, the neutrino-electron elastic scattering (ES) 
reactions at Super-Kamiokande and SNO, and the NC re- 
action at SNO could be affected. The SNO NC reaction 
is an axial current process , while the ES reaction de- 
pends on both axial and vector parts. Hence, the former 
is independent of the oscillation Hamiltonian while 
the latter is not. 

Since both and e^^ are strongly constrained, we 
set them to zero and vary ege, Cerj Ctt- Even with this 
reduction, the parameter space of the problem is quite 
large: different assignment of the diagonal and offdiago- 
nal NSI to electrons and u and d quarks yield different 
dependences of the oscillation Hamiltonian on the chemi- 
cal composition and different detection cross sections. To 
avoid complicating our main point with technical details, 
we limit our study to the case of NSI on quarks, assign- 
ing the same strength to the neutrino interactions with 
u and d quarks. 

For the solar neutrino analysis, we perform the stan- 
dard reduction of the 3x3 Hamiltonian to a 2 x 2 Hamil- 
tonian [12] . This involves performing a rotation in the 
/i — T subspace by the atmospheric angle ^23 and taking 
the first two columns/rows. This simplification is valid 
if (i) the 1 — 3 mixing angle is small: ^13 <C 1 and (ii) 
Gpne Cer^ ^™atm/(2^'!^)j with E^, being the neutrino 
energy and Am^^jj^ the difference of the squared masses, 
Am^jjj-j = to| — m\, as given by atmospheric neutrino 
data. The first requirement is ensured by the experi- 
mental bound from CHOOZ the second one can be 
checked to hold even for of order unity. The vacuum 
oscillation Hamiltonian then takes the usual form 



-A COS 26* A sin 20 
A sin 26* A cos 20 



(3) 



where A = Arn^ / [AE^) and Am? is the mass split- 
ting between the first and second neutrino mass states: 
Am? = m\ — m\. The matter contribution can be written 
(once again, up to an irrelevant overall constant) as: 

Gprie 



H, 



NSI 



V2 



1 + eii 

£12 



~12 



(4) 



where the quantities eij {i = 1, 2) depend on the original 
epsilons and on the rotation angle 023 ^ 



eil —fee — f'TT Sm 6^23 



£12 = -2 eer Siu 023 ■ (5) 



In Eq. (jSj), small corrections of order sin0i3 or higher 
have been neglected. We introduce a useful parameteri- 
zation: 



NSI 



A COS 2a ^e-^'"^ sin 2q; 



mat I Ae'^''t'sm2a -A cos 2a 



(6) 



Here the parameters A, a and (j) are defined as follows: 

tan2a = |ei2|/(l + en) , 2(f) = Argie^), 
A^Gpn, V[(l + eii)2 + ki2P]/2 . (7) 

In absence of NSI we have A = Gpnf.j \f2, a = and the 
Hamiltonian (jnj reduces to its standard form. 

Notice the appearance of the phase (f) in Eq. © . Since 
the phases of the basis states are chosen to make the 
elements in Eq. Q real, ^ cannot be simultaneously re- 
moved. This has been noted in the studies of terrestrial 
neutrino beams |0, , but overlooked in the solar neu- 
trino literature. 

What is the physical range for the parameters 0, a and 
(j>7 In the standard case a = 0, the physical range of is 
[0, 7r/2], including the so-called "light" and "dark" sides 
. A generalization to the NSI case is 

0e[O,7r/2], ae[0,7r/2], £ [-7r/2, 7r/2]. (8) 

These ranges of parameters cover all possibilities in 
Eqs. H3I6|) . For solar neutrinos, the range of (p could be 
cut in half, since points with (j) and —(f) give the same 
probability Pee that a Vg produced in the Sun is seen as 
a j/g in a detector. Moreover, the points (0, a,(/)) and 
(7r/2 — 0, 7r/2 — a, 0) are related by Pee ^ -P/i/i, which are 
equal in the 2-neutrino case by unitarity. 

Let us determine the expression for Pee. We first 
note that, because KamLAND selects Am^ > 10~^ eV^, 
coherence between the Hamiltonian eigenstates is com- 
pletely lost once one integrates over the neutrino energy 
spectrum and the neutrino production region in the Sun. 
The expression for the incoherent survival probability can 
be most easily derived in the basis that diagonalizes the 
matter Hamiltonian © . We obtain the familiar form 



Pee = [! + (!- 2Pc)cos20oCos20]/2, 



(9) 



where 9q is the mixing angle at the neutrino production 
point in the solar core and Pc is the level crossing proba- 
bility. These quantities contain all the effects of the NSI. 
The angle 00 is given by 



cos 20, 







cos 20 — xq cos 2a 
I + Xq + 2xq cos 20I.C 



(10) 



cos 20101 = sin 20 sin 2a cos 20 — cos 20 cos 2a. (11) 
Here xq = A/ A at the neutrino production point. 
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The expression for Pc is also easily found in the same 
basis, where it becomes apparent that the dynamics of 
conversion in matter depends only on the relative orienta- 
tion of the eigenstates of the vacuum and matter Hamil- 
tonians. This allows to directly apply the known analyt- 
ical solutions for Pc, and, upon rotating back, obtain a 
generalization of these results to the NSI case. For exam- 
ple, the answer for the infinite exponential profile Il8(.ll9l| 
A cx exp(— r/ro) becomes Pc = (exp[7(l — cos26'i.ci)/2] — 
l)/[exp{7) — 1], where 7 = AtttqA = nroAm'^ /E^. We 
further observe that since 7^1 the adiabaticity viola- 
tion occurs only when |0 — a| <^ 1 and_^~ which is 
the analogue of the small-angle MSW 0,1131 effect in the 
rotated basis. The "resonant" region in the Sun where 
level jumping can take place is narrow, defined by A ~ A 
|2l| . A neutrino produced at a lower density evolves adi- 
abatically, while a neutrino produced at a higher density 
may undergo level crossing. The probability Pc in the 
latter case is given to a very good accuracy by the for- 
mula for the linear profile, with an appropriate gradient 
taken along the neutrino trajectory. 



Pc ~ Q{A~A)e 



-7(cos2eroi-|-l)/2 



(12) 



where Q{x) is the step function, <d{x) = 1 for x > and 
&{x) — otherwise. We emphasize that our results differ 
from the similar ones given in in three important 

respects: (i) they are valid for all, not just small values of 
a (which is essential for our application), (ii) they include 
the angle 0, and (iii) the argument of the Q function does 
not contain cos 26*, as follows from "21]. We stress that 
for large values of a and (/i ~ 7r/2 adiabaticity is violated 
for large values of 6. 

Finally, to get an idea on the size of the day/night 
asymmetry, Aun = 2{N - D)/{N + D), (here D (N) de- 
notes the i^e flux at the detector during the day (night)) 
we can model the Earth as a sufficiently long (com- 
pared to the oscillation length) object of constant den- 
sity. For neutrino energies, this is appropriate for 
Am^ ^ 3 — 5 • 10^^ eV^. Introducing a small parameter 
X0 = A/ A, where A is evaluated for a typical density 
inside the Earth, we find, to the first order in xm, 



iDN 



sin 26'(cos 2a sin 26 + cos 2(/) sin 2a cos 26) 
' -[cos 26(0 (1 - 2Pc)]-^ - cos 26* 



(13) 

We verified that Eq. (|13|l gives a good agreement with 
precise numerical calculations for rig ~ 1.6 moles/cm"^. 
For the lower Am? region allowed by KamLAND, Am^ ^ 
1 — 3 • 10~^ eV^, the oscillation length is comparable to 
the size of the Earth, however, the averaging in Eq. H13|) 
still applies to a signal integrated over the zenith angle. 

In Fig. ^we plot the neutrino survival probability as a 
function of energy for several representative values of the 
NSI parameters. We take Am^ and 9 corresponding to 
the best-fit LMA point and choose the production point 
to be at r = O.li?0. Curve (1) is the standard interaction 
case, given for reference. The other three curves repre- 
sent the three qualitatively different regimes that are of 
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FIG. 1: The electron neutrino survival probability and the 
day/night asymmetry as a function of energy for Am^ — 7 x 
10~^ eV^, tan^ 9 = 0.4 and several representative values of 
the NSI parameters: (1) en = efi 
£11 = 4i = -0.008, £^2 = efa = -0.06; (3) e^i = 
-0.044, £5*2 = 42 = 0.14; (4) ej*! = efi = -0.044, 



42 = 0; (2) 

^11 — 



£12 



42 ~ —0.14. Recall that the parameters in Eq. JSJ equal 



interest to us. In the following we illustrate them in con- 
nection with observations. For definiteness, we consider 
real values of £12, both positive [(p — 0) and negative 
{(f) — 7r/2). As is clear from Eq. complex values 

(0 < 9!' < 7r/2) interpolate between these two cases. 



III. ANALYSIS OF DATA 

We now turn to the comparison of the NSI predictions 
with observations. To do this, we perform a best fit anal- 
ysis of the solar neutrino and KamLAND data along the 
lines ofRefs. |2^I3- In p articular, solar data include the 
radiochemical rates El Elllil , the SK ES zenith- 
spectra |2^ , the SNO day- night spectra 'so', "sT, mea- 
sured in phase-I and the SNO rates measured in phase- 
II 33]. For consistency, the NC rate prediction for SNO 
is treated as a free parameter because it is affected by 
an unknown change in the axial coupling of the quarks 
that could accompany the vector NSI considered in our 
analysis 34] . In our calculations, we use the updated 
BP04 35] Standard Solar Model (SSM) fluxes, electron 
density and neutrino production point distributions in 
the Sun. For KamLAND we considered the measured 
antineutrino spectrum with visible energies higher than 
2.6 MeV [33. 

The key ingredients of our analysis turn out to be the 
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rates and energy spectrum data from SNO and Super- 
Kamiokande. A comparison of the SNO CC rate with 
the Super-Kamiokande rate and the SSM indicates 
that, within the energy range accessible for the two ex- 
periments, the electron neutrino survival probability is 
about 30%. No other distinguishing features, such as a 
day/night asymmetry or spectral distortion, are seen at 
a statistically significant level js^ • In the case of the SM 
interaction, this turn out to be a very restrictive condi- 
tion; as seen in Fig. ^ the range of energies for which 
the survival probability is constant at 30% (henceforth, 
"the flat window") is barely large enough to cover the 
SNO energy window. On the low-energy end, the res- 
onant condition in the solar core increases the neutrino 
survival probability; on the high-energy end, the resonant 
condition in the earth causes a large D/N effect. Hence, 
values of the NSI parameters that "shrink" the flat win- 
dow, or shift it in the region disfavored by KamLAND, 
can be excluded. Conversely, if NSI increase the size of 
the flat window, new solutions may emerge. 



A. ei2 > 

A typical behavior for this case is exhibited by curve 
(3): the "step" in Pee becomes longer and the day/night 
asymmetry is not much smaller than in the SM case. 
These features point to a possible conflict with data. Our 
analysis confirms this expectation: a parameter scan Q 
for ei2 > reveals that a significant fraction of the pa- 
rameter space which is allowed by the accelerator-based 
data can be excluded by the solar/KamLAND data. 
As an example, we find that points with en = and 
> 0.14 (here and later, e^^ = e'^^ is assumed) are 
unacceptable at 90% confidence level (C.L.). If we keep 
A in the core of the Sun (at r — O.OSi?©) fixed to its 
standard value, A — Gprie/V^, we exclude points with 
^12 > 0-11 the same C.L. (for 1 degree of freedom, 
d.o.f, unless specified otherwise). The accelerator exper- 
iments allow values of order unity (in absolute value) for 
this parameter 0. We stress that the latter probe only 
|ei2|, while, as we show here, oscillation experiments are 
sensitive to the complex phase (p (or, for real epsilons, to 
the sign of eu). 



B. 



£12 < 



For ei2 close to zero (—0.08 < €"2 < 0), the only effect 
of the NSI is to flatten the part of the Pee curve around 
5-6 MeV, as illustrated by curve (2) in Fig. No new so- 
lutions appear and the allowed region in the plane 
is similar to that obtained with the SM interactions. This 
scenario has important implications for SNO, which can 
probe it by lowering its energy threshold. 

Finally, curve (4), obtained for e"i = —0.044, £"3 = 
—0.14, represents a novel and very interesting physical 
possibility. Its main feature is a significantly wider flat 
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FIG. 2; Regions of Am^ and tan^ 9 allowed at 
90, 95, 99, 99.73% C.L. (2 d.o.f) for SM interactions (left) and 
the NSI scenario (right) described by Eqs. 13161 . For the 
latter we used en = efi = -0.065, £12 = ^2 = -0.15. 



window, compared to the standard case. The key rea- 
son is the suppression of the day/night asymmetry on 
the high-energy end of the window. The physics of the 
suppression can be understood from Eq. H13|l . which, for 
(j) = 7r/2, gives Ajjn oc sm{29 — 2a). If the parameters 
are chosen in such a way that 9 and a in the Earth are 
comparable, the Earth regeneration effect is suppressed. 
Because of the difference in chemical composition, the 
difference 9 — a is larger in the Sun and, consequently, 
the evolution in the Sun is still adiabatic. 

A broader flat window allows the fit region to extend 
to lower values of Am^ compared to the standard case. 
While KamLAND excludes the middle part of the new 
region, the bottom part of that region (Am^ ~ (1 — 
2) X lO-'^ eV^) is, in fact, allowed H ^1- Thus, in 
addition to the usual two solutions, LMA-I and LMA-II, 
a completely new disconnected solution emerges, which 
we shall denote LMA-0. 

The situation is illustrated in Fig. El in which we 
compare the allowed regions in the standard case to 
those computed for chosen values of the NSI parame- 
ters: e^*! = -0.065, e5'2 = -0.15. The best-fit point 
in the LMA-0 region has Am^ = 1.5 x 10~^ eV'^ and 
tan^ e = 0.39, with = 81.7. For the same NSI pa- 
rameters, the has another minimum, — 79.9, at 
Am^ = 7.1 X 10~^ eV^ and tan^ 9 = 0.47, corresponding 
to the LMA-I solution. The quahty of the fit for LMA-0 
and LMA-I is comparable; if only the KamLAND rate, 
and not spectrum, information is used, the LMA-0 fit is 
slightly better (x^ = 73.0, against = 73.7 for the min- 
imum in the LMA-I region). For comparison, the best fit 
parameters for the standard case are Am^ = 7.1 x 10~^ 
eV^ and tan^ 9 = 0.43, with = 79.6 (using the Kam- 
LAND spectrum). 

We stress that the existence of the LMA-0 solution de- 
pends mainly on the value of a in the Earth and, to a 
lesser extent, on the value of the norm A. Hence, this 
solution persists for other choices of the NSI parameters, 
so long as they yield approximately the same a in the 
Earth. For example, if the diagonal interactions are as- 
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FIG. 3: The predicted KamLAND spectrum (top) and the 
time-averaged solar neutrino survival probability (bottom) 
for the LMA-0 best-fit point. For comparison, the standard 
LMA-I survival probability is also given. Refer to the text for 
details. 



sumed to be standard, e^i = 0, one finds a good fit in the 
LMA-0 region for e"^^ ~ -0.25. 

Our LMA-0 solution should not, of course, be confused 
with the "VERY-low-LMA" solution [s^, which arises 
under completely different physical assumptions, namely, 
if one assumes large (^ 5 — 8%) density fluctuations [Sg 
in the Sun [45l| . 

We note that the LMA-0 solution requires that the 
value of the (p angle be not too different from tt/2. Nu- 
merically, if we fix all the other parameters to the values 
of fig. 12 and vary ip, we find that LMA-0 disappears at 
90% C.L. for (f) < 0.457r. As is decreased further, the 
goodness of the overall fit decreases. Indeed, for = 
(positive £12), the survival probability has the features of 
curve (3) in Fig. which are disfavored, as discussed 
earlier. The LMA-I solution disappears at 90% C.L. for 
(j) < O.SItt. a scan over the region of £12 real and neg- 
ative gives exclusion of regions of the parameter space 
allowed by accelerator limits. For instance, points with 
£11 = and £"2 < —0.32 are unacceptable at 90% C.L. 
For A fixed to the standard value Gpne/V^ in the solar 
core, the limit is £5*2 < -0.19, at 90% C.L. 



Our choice of £"2 = —0.15 imphes £g^ = e^^ + £"^ ~ 
0.11 (see Eq. (O, in which we set 6*23 — 7i'/4), i.e., for ex- 
ample, £".,?' ~ £gT? ~ 0.05. This is about one order of mag- 
nitude smaller than the direct bound from CHARM j^. 
A more interesting question is whether the NSI param- 
eters of interest for the LMA-0 scenario could be tested 
with atmospheric neutrinos. For our specific case, the 
existing two-neutrino analyses |39l l40j | do not provide 
an answer, as the problem is essentially a three-fiavor 
one. Our investigation P| shows that regions exist in the 
space of the NSI parameters where the effect of NSI on 
the atmospheric neutrino observables is minimal and a 
satisfactory fit to the data is obtained. As an example, 
a point in this allowed region is £"g = £fg = —0.025, 
e^^ = ef^ = 0.11, e"^^ = ef^ = 0.08. 

The survival probabilities for the best-fit point of the 
LMA-0 solution and the standard LMA-I solution are il- 
lustrated in Fig. 13 (bottom). The curves represent proba- 
bilities averaged over time and over the production region 
inside the Sun for the *B and pep components of the solar 
neutrino spectrum according to |35j| . The probabilities 
for ^Be and pp neutrinos, not shown, are very close (with 
less than ~ 7% difference) to those for *B and pep respec- 
tively in the energy range of these neutrino fluxes. The 
energy intervals relevant to the different spectral compo- 
nents are also shown in the figure. The interval for ®B 
neutrinos is cut from below at E^, = 6.5 MeV; this ap- 
proximately corresponds to the threshold of T ~ 5 MeV 
in the electron energy at the SNO experiment. 

Interestingly, the LMA-0 solution has the features 
sought after in ^3|> where a sterile neutrino was intro- 
duced to eliminate the LMA-I upturn at SNO and im- 
prove the agreement with the Homestake rate. 

It is remarkable that, despite the wealth of data col- 
lected up to this point, such radically different scenar- 
ios as LMA-0 and LMA-I cannot be distinguished. The 
data expected in the next several years, on the other 
hand, should be able to resolve the ambiguity. First, if 
the SNO experiment lowers its energy threshold, it may 
be able, with sufficient statistics, to look for the upturn 
expected for the LMA-I solution. The absence of the up- 
turn would indicate the presence of NSI, or some other 
new physics. Second, the expected ''Be flux in the case of 
LMA-0 is lower, and the difference could be detected by 
the Borexino experiment (or by the future solar phase of 
KamLAND). Third, the small value of Am^ could be de- 
tected in the KamLAND spectrum data. The predicted 
spectra for LMA-I (standard interactions) and LMA-0 
are shown in Fig. O (top). It can be seen that the two 
are different at high energy where LMA-0 predicts more 
events. Thus, to make the discrimination it is necessary 
to both collect enough data and have a reliable calcula- 
tion of the antineutrino ffux for Ej^ >6 MeV. Finally, as 
evident from Fig. |3| the two solutions make dramatically 
different predictions for a pep experiment i46i| . 

While an observation consistent with the standard 
LMA-I solution would allow placing a very effective con- 
straint on the neutrino-matter interactions, a discovery 
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of a deviation consistent with the NSI signal would have 
truly profound particle physics implications. For exam- 
ple, according to Refs. [1,13' such interaction could be 
due to the operator of the form A1~'^Ir{H'^ aL){L<jH)lii 
(X v'^M^^{vv){lBl[i). For this operator to have an effect 
on the solar neutrino survival probability, the coefficient 
(X v^M'^ must not be too small, i.e., the scale of new 
physics M must not be much higher that the weak scale 
(Higgs vev v). Thus, by looking for the NSI signatures 
in solar/reactor neutrinos the experiments could in fact 
be probing new physics at the TeV scale. 

In summary, the present-day loose bounds on some of 
the neutrino interaction parameters introduce a serious 
uncertainty in the value of Ato^ extracted from solar and 



KamLAND data, allowing for a new, disconnected solu- 
tion. These uncertainties might be eliminated in the next 
several years, as more data are collected and analyzed by 
solar and KamLAND experiments. The constraints on 
the neutrino interactions presented here will be further 
extended. On the other hand, deviations from the SM 
neutrino interactions could indicate the presence of rad- 
ically new physics. We urge experimentalists to consider 
these points in their data analysis. 
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